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The large set of accurate data on differential cross section and analyzing power from the CERN 
LEAR experiment on pp — > n~n + in the range from 360 to 1550 MeV/c is well reproduced within 
a distorted wave approximation approach. The initial pp scattering wave functions originate from a 
recent AW model. The transition operator is obtained from a combination of the 3 Po and 3 Si quark- 
antiquark annihilation mechanisms. A good fit to the data, in particular the reproduction of the 
double dip structure observed in the analyzing powers, requires quark wave functions for proton, 
antiproton, and pions with radii slightly larger than the respective measured charge radii. This 
corresponds to an increase in range of the annihilation mechanisms and consequently the amplitudes 
for total angular momentum J — 2 and higher are much larger than in previous approaches. The 
final state 7T7t wave functions, parameterized in terms of 7T7t phase shifts and inelasticities, are also 
a very important ingredient for the fine tuning of the fit to the observables. 

PACS numbers: 12.39.Jh.; 13.75.Cs; 21.30.Fe; 25.43.+t 



I. INTRODUCTION 

The very accurate set of data from the LEAR experi- 
ment on pp — ► 7r~7r + measuring the differential cross 
section and analyzing power from 360 to 1550 MeV/c is 
still a challenge for theoretical models after more than a 
decade. Large variations are observed in the analyzing 
power Ao n as a function of angle at all energies, indicating 
presence of several partial waves already at low energies. 
However, recent model calculations 0, IE El IE IE lead 
to scattering amplitudes which are strongly dominated 
by total angular momentum J = and J = 1. The 
reason for this is the choice of a rather short range anni- 
hilation mechanism. The short range of the annihilation 
in the model calculations originates from the dynamics 
of baryon exchange in Refs. 0, IE El 0] or from required 
overlap of quark and antiquark wave functions for proton 
and antiproton in Refs. [E IE 13 • O n the other hand the 
experimental data on differential cross sections as well as 
those on asymmetries point to a significant J = 2, J = 3 
and even higher J contributions [EE! El ^2- 

All above mentioned models, for this reaction, use a 
distorted wave approximation (DWA). The ingredients 
for calculating the pp — > tt~tt + amplitudes consist of i) 
the initial pp scattering wave functions typ p (r) ii) a tran- 
sition operator 0(r',r) and iii) the final state 7T7t wave 
function $ n . ff (r / ). The complete scattering amplitude T 
itself, constructed in a DWA fashion, is 

T= Jdr'dr $^(r')0(r',r)% p (r). (1) 

For example in Ref. |5j] the transition operator 0(r',r) 
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was obtained from a combination of 3 Pg and 3 Si quark- 
antiquark annihilation model 

0(r', r) = N [Vs Po (r', r) + A Vs Sl (r', r) ], (2) 

where the relative strength A is a complex parameter and 
Nq an overall real normalization factor. In the same ref- 
erence ^E'pp(r) was provided by the 1982 Paris NN po- 
tential model and $ W7r (r') was a simple plane wave. 
This work did not succeed in reproducing the double-dip 
structure of the analyzing power and the forward peak 
in the differential cross section as seen experimentally 
at, for example, 497 MeV/c. All previously mentioned 
models exhibit similar difficulties. 

The aim of the present paper is to study possible im- 
provements of previous models. First of all, mesonic 
final-state interaction should be considered. The total 
energy of the pp — > 2-7T reaction for the studied data set 
is in the 2 GeV range. In this energy region the 7T7t in- 
teraction is characterized by several resonances ^4j. In 
Refs. 0, ||| the role of 7T7t final-state interactions was 
studied. In Ref. some improvement was obtained 
using a 7T7t model reproducing the real part of the 7T7t 
phase shifts with inelasticity parameters in all J ^ 
partial waves remaining close to 1. In Ref. |(| which ex- 
plores the 3 Pq part of the quark-antiquark dynamics in 
the transition operator, the final-state interaction, result- 
ing from meson exchange, affects mainly observables in 
the backward region. In both approaches the double-dip 
structure observed in experimental data of the analyzing 
power remains elusive and further study is still needed. 
Final-state interactions of two mesons in NN annihila- 
tion have also been studied, at quark level, within an 
extension of the quark rearrangement model @. Re- 
sults were reported for branching ratios of decays into 
various two-meson channels. Unfortunately there are no 
predictions from this work for differential cross sections 



2 



or analyzing powers. 

Within the approach of Ref. p| we will study the ef- 
fect of final-state interactions guided by the tttt coupled 
channel model of Ref. ^jj . Secondly we will study predic- 
tions following modification of the annihilation operator 
0(r', r). As remarked above this operator is rather short 
range in all present models. The reason for the short 
range of 0(r',r) is that in the quark- antiquark annihila- 
tion model the antiproton and proton have a relatively 
small radius since their quark wave functions describe 
only the qqq and qqq core ignoring the qq cloud. It could 
be a cause of discrepancy between theory and experi- 
ment. In Ref. it has already been noticed that an 
increase of the annihilation range improves substantially 
the theoretical description of the data. 

In the present work it is shown that the rnr final-state 
interaction is a very significant tool for the fine tuning 
of the fit to the observables. Furthermore the parame- 
ters that determine the size of protons and pions are also 
crucial. An increase of both proton and pion sizes, in 
closer agreement with their measured radii, allows for a 
much better fit to the experimental cross sections and an- 
alyzing powers. Expressions of the observables in term 
of the basic amplitudes together with the DWA ingre- 
dients are briefly recalled in Section II. The description 
of the final-state interaction is performed in Section III. 
The modifications of the range of the annihilation mech- 
anisms are studied in Section IV. Section V presents the 
final results and conclusions are summarized in Section 
VI. 



II. OBSERVABLES AND DWA INGREDIENTS 



and 



fj-1 



+ 1 tj+l 



xP'j(cos6). (8) 



The indices J and L = J ± 1 are the total and or- 
bital angular momentum of the pp system respectively. 
Pj(cos9) and Pj(cos8) denote a Legendre polynomial 
and its derivative. The angular momentum of the tttt sys- 
tem is = J. Because of parity conservation there are 
no L = J amplitudes in above expansion. Total isospin 
1 = for even J and 7 = 1 for odd J. In Eqs. and 
(|BJ p is the magnitude of the antiproton center of mass 
(CM) momentum. 

The partial wave amplitudes // for L = J ± 1 are 
calculated following the DWA method of Eq. One 
ingredient is the initial coupled spin-triplet ^ ^(r) wave 
function in configuration space as obtained in 17]. The 
operator 0(r,r') is constructed from the quark model 
description of protons and pions combined with the 3 Po 
and 3 Si quark-antiquark annihilation and rearrangement 
mechanism The last ingredient is the tttt scattering 
wave function ^^(r') which in this paper will be built ac- 
cording to a study of a realistic tttt scattering model ^| , 
while also comparisons will be made for a simple plane 
wave 7T7T final state. Subsequently one obtains the differ- 
ential cross section da/dVl and the analyzing power Ao n 
or left-right asymmetry for the proton target polarized 
normal to the scattering plane. 



The reaction pp — ^ 7r~7r + can be fully described in 
the helicity formalism by two independent helicity ampli- 
tudes F ++ (9) and F-\ (0). The angle 6 is the cm. angle 

between the outgoing ir~ and the incoming p. There are 
four possible observables ^(| 



\F + -\ 2 ), 



A 



On 



da 



Im(F ++ F*_), 



da 



(3) 
(4) 
(5) 
(6) 



So far only da/dQ and Ao n have been accurately mea- 
sured at LEAR l] . For completeness we recall here the 
partial wave expansion of the helicity amplitudes fl6| 



F ++W = ^Ev / (2^+l)/2 



Vj/y_! - VJTTfj +1 ] x Pj (cos 6), (7) 



III. FINAL-STATE INTERACTION 

First we assume no interaction between the final pions 
and describe the tttt scattering wave function $ W7r (r') as 
a plane wave. In this case the overall normalization of 
da/dil and the relative strength A between 3 Po and 3 Si 
annihilation mechanism arc parameters to be determined 
through x 2 minimization. Out of the twenty energies 
where polarization data are available jlj we chose a rep- 
resentative set of five energies, Ti a b = 66.7, 123.5, 219.9, 
499.2, and 803.1 MeV, corresponding to antiproton mo- 
menta of respectively pi a b = 360, 497, 679, 1089, and 
1467 MeV/c. Xi a b is the laboratory kinetic energy of the 
antiproton beam. At 123.5 MeV (pi ab = 497 MeV/c) 
and 219.9 MeV (p lab = 679 MeV/c) we do reproduce the 
results of using their parameters. Results for the set 
of five energies are shown in Figs. QHSl as short-dashed 
curves. The overall fits are poor with exception of the 
analyzing power A 0n at 803.1 MeV (p\ ah = 1467 MeV/c). 
The double-dip structure of Aq u is not reproduced at the 
three lowest energies of the set considered here, which 
can be attributed to the very small values of the J > 2 
amplitudes predicted by this model. A lack of substan- 
tial J > 2 amplitudes is also evident in the predictions 
for da / d^l at all energies. The forward peak is poorly 
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FIG. 1: Differential cross section and analyzing power of the reaction pp — > n n + at Ti a b = 66.7 MeV (pi a b = 360 MeV/c) 
Experimental data from Hasan et al. dj. The different curves are described in the text. 



200 



150 



: s 

CO 



3 

c 



100 



50 




1.0 -0.5 0.0 0.5 

cos(0) 



1.0 

7, 




FIG. 2: As in Fig.[T]but for T lab = 123.5 MeV (pi ab = 497 MeV/c) 



reproduced, and also the backward peak, prominent in 
the data at Ti a b = 499.2 MeV, is missing in this simpli- 
fied model. Similar findings were obtained previously at 
Ti a b = 123.5 and 219.9 MeV by [gj. As shown in the 
analysis of ^l| we recall that large J > 2 amplitudes are 
needed to explain the angular dependence of Aq u and 
da / d£l . 

One possibility to enhance the amplitudes of higher J 
values could be to introduce irir final-state interaction. 
The elastic 7T7t — > 7T7t amplitude is known from threshold 



up to the total relativistic 7T7t energy yfs = 1800 MeV 
mainly from analysis of the irN — > ttttA?' reaction. The 
extracted 7T7t — > 7T7t amplitudes can be parameterized in 
terms of phase shifts Sj and inelasticities r\j where J = 
0, 1, 2 and 3 0. In a coupled channel model of 
7T7r, KK and pp was proposed to reproduce these phase 
parameters Sj and r] j for J = 0, 1, 2 and 3. The $ WT (r') 
wave functions required for the final-state interaction of 
pp — > 7T7T, can be constructed from this coupled channel 
model in a straightforward manner. However, since the 
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FIG. 3: As in Fig.[T]but for T lab = 219.9 MeV (pi ab = 679 MeV/c). 




FIG. 4: As in Fig.[]]but for T lah = 499.2 MeV (pi ab = 1089 MeV/c). 



needed energy range in pp — > tttt is from y/s = 1910 to 
2272 MeV one has to rely on extrapolation of the coupled 
channel model results beyond y/s — 1800 MeV. Calcula- 
tions with the corresponding ^^(r') show a significant 
sensitivity of the observables to the tttt final- state in- 
teraction. Even so the tttt scattering amplitude from the 
extrapolated coupled channel tttt model still does not im- 
prove the predictions of pp — > tttt obtained previously 
with 7T7T plane waves. But in this study we did find that 
the off-shell part of $ W7r (r') does play a very minor role. 
In particular it was observed that the same predictions of 



observables da/dQ and Ao n can be obtained using only 
the asymptotic part of the tttt wave functions. We then 
exploit this fact by asserting for the remainder of this 
paper that in pp — > tttt the unknown tttt final-state scat- 
tering can be fully described by the asymptotic tttt wave 
functions parameterized by the tttt phases Sj and inelas- 
ticities r\j. 

We then avoid the problem of extrapolation to high 
y/s by including these new parameters Sj and rjj in the 
minimization process to obtain realistic fits to da/dfl and 
Ao n . Results of this second fitting procedure are shown 
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FIG. 5: As in Fig. □ but for T la . h = 803.1 MeV (pi ab = 1467 MeV/c) 
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FIG. 6: Various fits for Ti a b = 219.9 MeV (piab ~ 679 MeV/c). The solid curve is as in Fig. El The long-dashed curve represents 
a fit with final-state interaction, where a = 1.09 fm -2 and (3 = 1.51 fm -2 . The short-dashed curve is a fit with no final-state 
interaction, where a — 1.27 fm -2 and /3 — 1.53 fm -2 . 



as the long-dashed curves in Figs.mHl Switching on final 
state 7T7r interaction improves the fit of A$ n by readjust- 
ing the strength of the helicity amplitudes of different J. 
This allows for prediction of a double-dip structure at 
lower energies. This feature is a crucial requirement of 
the data and shows the need for incorporating the final- 
state interaction of the pions in some form. However 
the predictions for da/dQ show only a modest improve- 
ment over the model without final-state interaction and 



the question arises whether there is additional freedom 
within the model to ameliorate the present fit. So far, 
the only variable parameter in the annihilation operator 
is the relative strength A unless one allows variations of 
the range of the annihilation mechanism controlled by 
the parameters a and (3 Q. In the next section we will 
investigate the effects of variations in a and (3. 
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IV. MODIFICATION OF THE ANNIHILATION 
RANGE 

In order to derive the annihilation operator 0(r',r) 
@, one can describe the proton, antiproton and pions 
in terms of quarks and antiquarks with the use of Gaus- 
sian wave functions Q. This amounts to approximate 
quark confinement by solving the Dirac equation with 
either a scalar or vector harmonic oscillator potential. 
The proton (antiproton) intrinsic wave function for the 
annihilation mechanism is : 



tp N {Yi,r 2 ,r 3 ) = N N exp 



xxjv (spin, isospin, color), (9) 

where are the quark (antiquark) coordinates and r N 
the nucleon (antinucleon) coordinate. An S-wave meson 
intrinsic wave function is: 



,(ri,r 4 ) = N M exp 



ft 



i=l,4 

X% M (spin, isospin, color). (10) 

Here ri and r 4 are the quark and antiquark coordinates, 
respectively. The coordinate of the meson is r M . 

Typical parameter values used before are a = 
2.80 fm~ 2 and ft = 3.23 fm -2 which correspond to a pro- 
ton (antiproton) radius of 0.60 fm and a meson radius of 
0.48 fm [19j. This value of a describes the qqq (qqq) core 
of the proton (antiproton) while the measured charge ra- 
dius, which for the proton is about 0.8 fm, includes also 
the mesonic cloud. Explicit expressions in terms of a and 
ft for the transition operators V3p (r',r) and V~3 Sl (r',r) 
of Eq. (|2J) can be found in 5j. However, one can also 
argue that in modeling the Gaussian wave functions as 
in Eqs. © and l|10|) . values for a and ft simply should 
be in accordance with known charge radii of the proton 
and pion. The measured pion charge distribution radius 
[13 is (r 2 ) 1/2 = 0.663 ± 0.006 fm. For the proton we 
find (r 2 ) 1 / 2 = 0.870 ± 0.008 fm in the literature 03 . We 
are thus left with a certain freedom when it comes to 
choosing values for a and ft and we can wonder about 
the effects on the annihilation mechanism. The values of 
the parameters a and ft determine effectively the range of 
the annihilation mechanism. Let us assume we increase 
the size of the proton and the pion from their original 
values in [j| by decreasing a and ft. Then the integration 
over the intrinsic quark coordinates of, for example, the 
wave functions Vp(i"i, r 2 , r 3 ), ip p (r 4 , r 5 , r 6 ), ^+(ri,r 4 ), 
and 4>n- (i"2, 1*5) will result in a larger quark overlap. The 
corresponding annihilation operator 0(r',r) will have a 
longer range and therefore higher partial waves will con- 
tribute to the total amplitude T in Eq. (Q) as re q uired b y 
the analyses of experimental data of Refs. [9l flflt fill Il2l | . 

Changes of a and ft can furthermore be linked to rela- 
tivistic corrections of the pion wave function </>M(i"i,r 4 ) 



as the outgoing energy of the pion is much larger than 
its rest mass. In practical terms this means that in the 
center of mass, where the calculation is performed, the 
pion wave function should not be described as a sphere 
anymore. A proper treatment requires a Lorentz boost 
of the pion intrinsic wave function from its rest frame to 
the CM frame. Thus one expects a change in the geom- 
etry of the Gaussian shape of the pion |2l| . This affects 
the overlap integral which could be mocked by a si- 
multaneous alteration of a and ft. We therefore take a 
and ft to be variable parameters to the extent that they 
still satisfy physical conditions. 



V. FINAL FIT 

The results of our final fit which includes now values 
of the parameters a, ft, \ — | A| exp(z6>A), 5j, rjj, and 
of the overall normalization Nq are shown as solid curves 
in Figs. The improvement over the previous fits 

is dramatic but requires both increased sizes of the an- 
tiproton, proton and pions as well as a tuned final-state 
interaction. The main achievement is the reproduction 
in the differential cross section da/dVl of the character- 
istic forward peaks at lower energies (Ty, = 66.7, 123.5 
and 219.9 MeV) and backward peaks at higher energies 
(TLb = 219.9 and 499.2 MeV). This clearly indicates that 
the increase of the annihilation range now produces much 
larger amplitudes for J = 2 and higher. The predictions 
for the double-dip structure of the analyzing power Ao n 
compare much better with the experimental results espe- 
cially for T lab > 219.9 MeV. 

In addition to the double-dip structure of Ao n the ex- 
perimental data display another characteristic feature: 
the asymmetry shifts from predominantly negative val- 
ues at lower energies toward positive values at higher 
energies. Our final fit accounts for this pattern. The 
quark model parameters a, ft, A, and the overall norm 
No resulting from this fit are listed in Table |U and the 
phase shifts Sj and inelasticities 77,7 of the final-state in- 
teraction with their dependence on y/s in Table ITT1 Note 
that the latter have been readjusted in the final fit and 
differ from the Sj and r/j obtained in the fit presented 
in section UTTl The irir phases for J = 0, 1 are small but 
for J = 2, 3, 4 they are substantial. The inelasticities 
r/j, in particular for J = 2 may indicate the presence of 
resonances in this channel. From Table [I] it is also clear 
that the parameters a and ft are almost constant with 
energy but have lower values at 219.9 MeV. Neverthe- 
less, in Fig. we show (long-dashed line) that one can 
obtain a fit of similar quality keeping the size parame- 
ters a = 1.09 fm~ 2 and ft = 1.51 fm~ 2 close to the ones 
at the other four energies. The relative strength A ex- 
hibits an energy dependence. This dependence shows a 
smooth decrease of |A| with increasing y/s, which indi- 
cates that the 3 Pq mechanism preponderates at higher 
energies. The phase 0\ of A should be compared with 
the value 9\ = 180° of [5j. Furthermore, we find that 
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FIG. 7: Predictions at TLb = 66.7 MeV (pi a b = 360 MeV/c) for the observables Ai a and ^4 SS obtained with the final fit. 



a decreases by about a factor 2.3 while f3 decreases by 
about 2.2. In other words, the proton size in this model 
increases to (r 2 ) 1 / 2 = 0.91 fm and the pion radius is now 
(r 2 .) 1 / 2 = 0.71 fm which is within 7% of the values of ref- 
erences 0| and [2(j mentioned in the previous section. 

The dramatic improvement that occurs when the size 
parameters a and (3 take on smaller values, tends to mask 
the equally important role of the final-state interaction. 
If the final-state interaction is turned off from the very 
beginning by fixing the phase parameters Sj = and 
r/j = 1 while a, (3, A are allowed to vary, a and f3 again 
decrease significantly, which confirms an increased anni- 
hilation range within our model. Nevertheless, the result- 
ing fit without final-state interaction is far from satisfac- 
tory and it is only when we include final-state interaction 
that we recover the fit quality discussed earlier in this sec- 
tion. This is illustrated at T lah = 219.9 MeV (p lab = 679 
MeV/c) in Fig. EI Tbe short-dashed curve represents the 
fit without final-state interaction, and for which a = 1.27 
fm -2 and (3 — 1.53 fm~ 2 . This short-dashed curve in 



Ti ab [MeV] 


66.7 


123.5 


219.9 


499.2 


803.1 


piab [MeV/c] 


369 


497 


679 


1089 


1467 


a [fm" 2 ] 


1.20 


1.19 


1.11 


1.20 


1.20 


P [fm- 2 ] 


1.54 


1.51 


1.00 


1.54 


1.54 


\M 


1.200 


1.064 


0.662 


0.545 


0.424 




197.10 


125.68 


166.83 


176.32 


183.29 


N [10 s ] 


5.251 


4.804 


5.046 


4.222 


4.667 



TABLE I: Quark model parameters as function of Ti a b . 



Fig. H3 should be compared with the short-dashed curve 
in Fig.Elfor which a = 2.80 fm" 2 and (3 = 3.23 far 2 . 
The differential cross section has improved significantly 
in the forward hemisphere. However the analyzing power 
is only marginally better. The solid line of Fig. is with 
inclusion of the final-state interaction. Similar results 
are obtained at the other energies, which confirms that 
in order to reproduce the LEAR data both an increased 
annihilation range as well as interaction in the final 7T7t 
state is necessary. 

With the present parameters, we can proceed to com- 
pute other spin observables for the pp — ► tt~it + reaction. 
We take the opportunity to present for each energy con- 
sidered above the predictions for the spin observables Ag s 



Tiab [MeV] 


66.7 


123.5 


219.9 


499.2 


803.1 


[MeV] 


1910. 


1937. 


1983. 


2111. 


2242. 


Vo 


0.840 


0.989 


0.937 


1.00 


1.00 


So 


-4.57 


-4.37 


-6.59 


-5.24 


-0.63 


m 


0.998 


0.995 


0.921 


0.999 


0.914 


Si 


8.05 


4.89 


-4.14 


-3.24 


-1.88 


V2 


0.756 


0.669 


0.760 


0.818 


0.309 


52 


35.74 


24.49 


55.57 


37.66 


74.95 


V3 


1.00 


1.00 


0.808 


1.00 


1.00 


8 3 


59.66 


41.41 


62.45 


48.15 


-19.59 


m 


1.00 


0.997 


0.858 


1.00 


0.990 


5i 


36.93 


-41.29 


6.19 


-60.70 


61.77 



TABLE II: Phases shifts 5j and inelasticities r\j of the final- 
state interaction for < J < 4 as function of y/s. 
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FIG. 8: As in Fig.EJbut for T lab = 123.5 MeV (pi ab = 497 MeV/c 

1.0 





0.5 



jr o.o - 



-0.5 



-1.0 




-1.0 - 



0.5 0.0 
cos(0) 



0.5 1.0 




-0.5 0.0 0.5 
cos(0) 



1.0 



T, ^— 219.9 MeV 

lab 



FIG. 9: As in Fig.|71but for T lab = 219.9 MeV (pi ab = 679 MeV/c). 



and A ss introduced in Eqs. JSJ) and We remind the 
reader that the spin observables are related by 

^0n + -A$ s + A 2 SS = 1 , (11) 

and that the reaction pp — > ir~ir + has only three indepen- 
dent observables within a sign ambiguity, and of course 
da/dil is chosen as one of them. At each energy the cor- 
responding parameters of Tables HI and ITU are used. The 



results are shown in Figs.l7HrTl The spin observables A? s 
and A ss exhibit again the typical structure with two ex- 
trema as a function of angle, which take often the form 
of a double dip. There is as yet no data to compare with. 
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FIG. 10: As in Fig.EJbut for Ti ab = 499.2 MeV (pi ab = 1089 MeV/c 
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FIG. 11: As in Fig.^but for T lab = 803.1 MeV (pi ab = 1467 MeV/c). 



VI. CONCLUSION 

The extensive set of data of differential cross sections 
da/dfl and analyzing powers Aq u from the LEAR 
experiment on pp — > 7t~tt + in the range pi ab = 360 to 
1550 MeV/c can be fitted by the combined mechanisms 
3 Po and 3 Si of the quark-antiquark annihilation model. 



The initial pp relative wave functions were taken from 
[l?| . It is important to include the final state nn in- 
teraction and employ quark wave functions for proton, 
antiproton, and pions with radii which are slightly larger 
than the respective measured charge radii. Previously 
used hadron intrinsic quark wave functions Q describe 
only the quark core of the hadrons without the qq cloud 
and their parameters therefore correspond to a consid- 
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erably smaller radius. Increased hadronic radii lead to 
an increase in range of the annihilation mechanism and 
as a result amplitudes for J = 2 and higher are much 
larger than before. This feature in the model is essen- 
tial since the experimental data on A 0n exhibit a double 
dip, which indicates the presence of substantial ampli- 
tudes of J = 2 and higher, already at lower momenta 
piab • The relative strength of the 3 P and 3 Si mecha- 
nisms shows a smooth energy dependence and suggests 
that the 3 Po mechanism becomes more dominant at the 
higher energies. It is however noted that the pronounced 
forward and backward peaks in the cross section require 
the presence of both 3 P and 3 Si mechanisms. 
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